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We study the S-wave interaction of mesons with baryons in the strangeness S = sector 
in a coupled channel unitary approach. The basic dynamics is drawn from the lowest order 
pr) | meson baryon chiral Lagrangians. Small modifications inspired by models with explicit 

vector meson exchange in the t— channel are also considered. In addition the ttttN channel 
is included and shown to have an important repercussion in the results, particularly in the 
I = 3/2 sector. The A* (1535) resonance is dynamically generated and appears as a pole in 
the second Riemann sheet with its mass, width and branching ratios in fair agreement with 
experiment. A A(1620) resonance also appears as a pole at the right position although with 
a very large width, coming essentially from the coupling to the ttttN channel, in qualitative 
agreement with experiment. 
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1 Introduction 



The introduction of effective chiral Lagrangians to account for the basic symmetries of QCD and its 
application through xPT to the study of meson meson interaction or meson baryon interaction 
0, |3|, (|, Q has brought new light into these problems and allowed a systematic approach. Yet, 
xPT is constrained to the low energy region, where it has had a remarkable success, but makes 
unaffordable the study of the intermediate energy region where resonances appear. In recent 
years, however, the combination of the information of the chiral Lagrangians, together with the 
use of nonperturbative schemes, have allowed one to make prediction beyond those of the chiral 
perturbation expansion. The main idea that has allowed the extension of xPT to higher energies 
is the inclusion of unitarity in coupled channels. Within the framework of chiral dynamics, the 
combination of unitarity in coupled channels together with a reordering of the chiral expansion, 
provides a faster convergence and a larger convergence radius of a new chiral expansion, such 
that the lowest energy resonances are generated within those schemes. A pioneering work along 
this direction was made in || |7|, |8j where the Lippmann Schwinger equation in coupled channels 
was used to deal with the meson baryon interaction in the region of the A?"* (1535) and A(1405) 
resonances. Similar lines, using the Bethe Salpeter equation in the meson meson interaction, were 
followed in || and a more elaborated framework was subsequently developed using the Inverse 
amplitude method (IAM) [|Hj] and the N/D method [jlTJ. The IAM method was also extended 
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Figure 1: Diagram representation of the Bethe-Salpeter equation. 



to the case of the meson baryon interaction in and in ]T3[ where a good reproduction of the 
A resonance was obtained using second order parameters of natural size. The N/D method has 
also been used for ttN scattering []rjj and in the K~N and coupled channels system in A 
review of these unitary methods can be seen in |16[]. The consideration of coupled channels to 



study meson baryon interactions at intermediate energies has also been exploited in |T7| using the 
K-matrix approach, although not within a chiral context. 

The Bethe Salpeter equation was also used in the study of the the meson baryon interaction in 



the strangeness S = — 1 sector in |18| and in the S = sector around the iV*(1535) region in |19 



In this latter work, aimed at determining the N*N*n coupling, only the vicinity of the resonance 
was studied and no particular attention was given to the region of lower energies. Subsequently a 
work along the same lines using the Bethe Salpeter equation, but considering all the freedom of 
the chiral constraints, was done in [20] and a good reproduction of the experimental observables 
was obtained for the 1 = 1/2 sector. Those works considered only states of meson baryon in the 
coupled channels and both in [|T^] and |2(| the ttttN channel was omitted. This channel plays a 
moderate role in the 1=1/2 sector |21| but, as we shall see, it plays a crucial role in the I = 3/2 
sector. 

Our aim in the present work is to extend the chiral unitary approach to account for the ttttN 
channel, including simultaneously some other corrections inspired by vector meson dominance 
(VMD) which finally allow one to have a reasonable description of the meson baryon interaction up 
to meson baryon energies of around 1600 MeV. The iV*(1535) resonance is generated dynamically 
in this approach and the mass, width and branching ratios are obtained in fair agreement with the 
experiment. The phase shifts and inelasticities for irN scattering in that region are also evaluated 
and good agreement with experiment is also found both in the 1=1/2 and 1 = 3/2 sectors. In 
addition some trace of the A (1620) is found, linked to the introduction of the 7T7riV channel, with 
a pole in the second Riemann sheet with the right energy, albeit a large width. 

A side effect of the calculations is that we determine the s-wave part of the 7riV — > nnN 
transition amplitude, revising previous determinations in [23] and [24], and together with the P- 
wave amplitudes previously determined, we obtain a good reproduction of the cross sections for 
these reactions. 



2 ttN scattering in a 2-body coupled channel model 
2.1 Basic 2-body model 

In this section we study the meson baryon scattering in S-wave in the strangeness S = sector. 
We shall make use of the Bethe Salpeter equation in coupled channels considering states of a 
meson of the 0~ octet and a baryon of the l/2 + octet, as required by the SU(3) chiral formalism. 
For total zero charge we have six channels, ir~p, 7r°n, rjn, K + Y>~ , i^°S°, and K°A. 
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Table 1: Cy coefficient in the potential. = Cji 



K + E~ if u E u K U A n-p ?r u n rjn 

if A - x - :; - ;; 



2 2 2 

1 -V^ 



7T°n 

rjn 



The Bethe Salpeter equation is given by 



T = V + VGT 



(1) 



where G is the product of the meson and baryon propagators. The diagrammatic expression of 
the T matrix is shown in Fig.|l|. Following we take the kernel (potential) of the Bethe Salpeter 
equation from the lowest order chiral Lagrangian involving mesons and baryons 



MB^M'B' 



Tr 



(2) 



where B and are the SU(3) matrices for the octet baryon field and the octet meson field 
respectively, and / is the weak decay constant of the meson. From this Lagrangian, the transition 
potentials between our six channels are given by 



(3) 



with the initial (final) baryon spinor u(p) (u(p')) and the initial(final) meson momentum k (k r ). 
The coefficients Cy, reflecting the SU(3) symmetry of the problem, are obtained from eq. (0) and 
shown in Table 0. 

We are interested in the value of the T matrix element for the on-shell meson-baryon system 
at a certain center of mass energy y/s or P = p + k = (y/s, 0, 0, 0), and we express it as T(y / i). 
Yet, in the diagrammatic expression of the Bethe Salpeter equation, Fig.Q, one can see that the 
V and T matrices in the second diagram on the right hand side, can be (half) off-shell since the 



BS equation is an integral equation. However, it was shown in [18| that the off shell part of these 
matrices in the loops could be incorporated into a renormalization of the lowest order Lagrangian 
and hence only the on shell parts are needed. Thus, they factorize out of the integral and the 
original BS integral equation is then reduced to an algebraic equation which is easily solved with 
the solution in matrix form 



T(V~s) = [V( v rs)- 1 -G( v rs)]- 1 or T(v^)^ 1 = V(v^)- 1 - G(v^) 



(4) 



where the G function is a diagonal matrix representing the loop integral of a meson and a baryon. 
The i-th element is thus expressed as 



Gi(P) 



d 4 q 



2Mt 



1 



(2tt) 4 (P - q) 2 -Mf + ieq 2 - mf + ie 



(5) 
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with Mi the baryon mass and rrii the meson mass. 



In the present (Mandl and Shaw) normalization [p2|, the T matrix is related to the S matrix 



as 



Syft, P 3 ) = 1 - i (2n)W(P l - 1 l Ml 



V 2 y Ei\ Ej V 2ui y 2^- 



Tij(Pii Pj) 



(6) 



where V is the volume of the normalization box, Ei (Ej) and Ui (tUj) are the energies of the 
incoming (outgoing) baryon and meson respectively. In this notation, the unitarity condition for 
any partial wave amplitude (in particular the S-wave which is our case) is given by 



where 



J(s - (Mi + mi) 2 ) (s - (Mi - mi) 2 ) 
Q.(v^) = ^ — ^= — —0(V^ - (Mi + rm)) 

is the on shell center of mass momentum of z-th meson-baryon system. 
The matrix elements of the potential which we substitute into eq. M) is 



(7) 



(8) 



fy(VS) = -C llT ^(2V7s -M<- Mj )J m + E ^ S) 



\ 



2M< 



(9) 



where we introduce different weak decay constants for different mesons. We use the values 

f n = 93 MeV , f K = 1.22U , /„ = L3/ w (10) 

taken from %PT [[[[]. 

The meson-baryon 2-body propagator which we substitute into equation M) is 



„ , rs 2Mj f / \ i mf Mf — m 2 + s M 2 



+ 



(4tt) 



2s 



log (s - (M? - m^) + 2v^Qi(v^)) + log (s + (M? - m, 2 ) + 2^~sQi(V~s 



log (-3 + (M 2 - m t 2 ) + 2yTsQiWa)) - log (-a - (M, 



2 2\ 

, — m? 



2^Q t (V^)) 



which is the 1-loop integral fl5|) done with dimensional regularization. Its infinity is canceled 
by higher order counter-terms. The first terms cij(/i) are real constants and stand for the finite 
contribution of such counter-terms. We treat these a«(/i) as unknown parameters and determine 
them from fits to the data. The imaginary part of the above propagator is 



ImiG^)} 



MjQijJs) 



(12) 



which leads, using eq. (J|), to eq. (0). Therefore, in the present model, unitary is exactly fulfilled. 

In order to keep isospin symmetry in the case that the masses of the particles in the same 
multiplet are equal, we choose cii(fi) to be the same for the states belonging to the same isospin 
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Figure 2: Propagators of the meson-baryon systems using the parameters of eq flT3"|). The K°TP 
and TT°n propagators are omitted because they are almost the same as those of K + T,~ and n~p 
respectively. The continuous(dashed) lines stand for the real(imaginary) part. 



multiplet. Hence we have four subtraction constants a^-A^/i), a^^), aK\{^) and a^s(/x). A best 
fit to the data with eq. leads us to the following values of these parameters 

IM = 1200 MeV, a^Gu) = 2.0, a vN {fi) = 0.2, a KA {fi) = 1.6, a X s(/i) = -2.8 . (13) 

In Fig.0, we show the G matrix elements which are obtained with this set of parameters. These 



propagators are essentially the same as those obtained in |1| where the regularization was done 
using a three- momentum cutoff (\q\ < 1 GeV) 



Re[G i { v G)]=V 



IGeVtfq Mi 



1 



(2tt)3 E(q) yfi - E(q) - u(q) 2u{q) 



(14) 



with E(q) 



M? + q 2 , u(q) = \Jm 2 + q 2 and 

a nN = 35 MeV, a vN = 16 MeV, a KA = 40 MeV, ~a KYl 



-21 MeV 



(15) 



where V means the Cauchy principal value. 

The resulting T matrix elements of Sn (isospin 1/2) and S31 (isospin 3/2) tcN elastic scattering 
are shown in Fig.[| In these graphs we plot the quantity 



MjQ^s) M 3 Q 3 {JS) 



Tij(Vs) 



(16) 



in order to compare with the data of the CNS analysis [P5[| . We find a qualitative agreement 
with the data in the energy range from threshold to 1600 MeV. In the figure of the experimental 
analysis one can see the manifestation of the resonances A r *(1535), iV*(1650) in the Sn amplitude 
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Figure 3: Scattering amplitude for the Sn and S 3 i irN partial waves. The left and right figures 
stand for the amplitudes calculated in our model and the data analysis |25| respectively. The 
continuous(dashed) lines stand for the real(imaginary) part. 



and A(1620) in the Ssi one. The calculated Sn amplitude also exhibits a resonance structure 
around 1535 MeV. The generation of this resonance is common to all the unitary chiral approaches 
H [19|, |20|. In section 4, after we include new elements in the theory, we shall investigate this 
resonance by searching for poles in the second Riemann sheet of the complex plane. 

At energies beyond 1600 MeV, the calculated amplitudes are qualitatively different from the 
data, and the iV*(1650) and A(1620) do not show up. In the philosophy that the resonances 
obtained using the lowest order Lagrangian and the present unitary scheme are simply meson 
baryon scattering resonances, (qualifying for quasibound meson baryon states [16|]), the nongener- 
ation of a particular resonance would indicate that it is mostly a genuine state (approximately a 
3q system). However, such resonances could be also obtained in a unitary approach provided one 
used information related to this resonance which would be incorporated in the higher order La- 
grangians. For the case of the meson meson interaction it is known |26| that the 0(p A ) Lagrangian 
is saturated by the exchange of vector meson resonances, which are not generated in the BS ap- 
proach of |§. Actually, in J2(| the iV*(1650) resonance is also reproduced by introducing counter 
terms which effectively account for higher order corrections, much in the way as the (genuine) 



p resonance was reproduced in the study of the meson meson scattering in |27| , p8| . As quoted 
before, the agreement below 1600 MeV is only qualitative. Indeed, both the real and imaginary 
parts of the Sn amplitude are somewhat overestimated in the theory, and so is the case for the 

amplitude where the theoretical imaginary part clearly overestimates the experimental one. 

The phase-shifts and elasticities rji(yfs) are given by 



1 

- arccos 
2 



Re 



1 - 2i 



T- 

1 it 



x sign [Re [-Tu 



1 - 2i—^T u 



(17) 
(18) 
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Figure 4: Phase-shifts and inelasticities for Su and S31 nN scattering. The continuous(dotted) 
lines correspond the calculations(data analysis p5|). 



Fig.|4] shows the phase-shifts and inelasticities. In the phase-shifts graph we can see again the 
qualitative agreement at energies below 1600 MeV. On the other hand, the inelasticities are not 
reproduced below the first open meson baryon threshold. In this 2-body model the threshold of 
inelastic scattering for the Su case is the rjn threshold which appears at 1487 MeV. Below this 
energy, the calculated inelasticities are zero and do not agree with the data. This situation is 
even clearer in the S31 case. The threshold of inelastic scattering is in this case the KH threshold 
which appears at 1690 MeV. The big inelasticities observed in the data below that energy are not 
reproduced in the present 2-body approach. The only inelastic channel opened below that energy 
is the 7T7riV channel, and the experimental data is telling us that the influence of this channel in 
the S31 amplitude should be very important. We will include the 7nrN channel in the next section. 



2.2 Improved 2-body model 



The vector meson dominance (VMD) hypothesis is phenomenologically very successful. In this 
hypothesis, the meson baryon interaction is provided by vector meson exchange. For example, 
the tiN — > ttN process is described by p meson exchange in the t-channel, as shown in Fig.|5]. 
It is interesting to note that the result from p exchange provides the amplitude obtained from 
the lowest order chiral Lagrangian. For example, in the case of 7i~p — > n^p, p exchange in the 
t— channel gives (see 



29] for the pNN coupling within the VMD hypothesis) 



.m v G v 



q 2 — ml 



.m v G v 



\k + k') 



—1 



4/ 



(19) 
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Figure 5: The mechanism of — > ttN interaction in the vector meson dominance hypothesis. 
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Figure 6: Correction factor of the chiral coefficient. 



which follows from the relation G\j ' f 2 = 1/2, which is a result of the VMD hypothesis. This 
indicates that the lowest order chiral Lagrangian is an effective manifestation of the VMD mech- 
anism in the vector field representation of the vector mesons. According to this consideration, we 
introduce in the chiral coefficient a correction to account for the dependence on the momentum 
transfer of the p propagator. Thus we replace 



C; 



Cij x 



dk' 



-m: 



4tt (k'-k) 2 



at yG > 



(20) 



where ysf^ is the energy where the integral of eq. fl2HD is unity, and which appears in between 
the thresholds of the two i,j channels. At very low energies where xPT is used, this correction 
is negligible but this is not the case at the intermediate energies studied here. For example, the 
correction for n~p — > n~p element is calculated as 



dk' 



—m r 



Atx (k' - k)< 



m 2 m 2 + 2k°k'° + 2kk' - m 
rylog : 



mz 



Akk' 



m 2 p + 2k°k'° 



2kk' 



mi 



mi 



(21) 



with m p = 770 [MeV], and is shown in Fig.||| left. One can see that p meson tail reduces the 
coefficient about 25% at energies around 1500 MeV. Similarly, for the strangeness exchanging 
process, we consider K* exchange in the t-channel. For example, the correction of 7r - p — > K°A is 
calculated with mx* = 892 MeV and shown in Fig.[| right. 

With this correction, after retuning the subtraction constants aj(/i), we obtain the T matrix 
shown in Fig.|7|. As we can see, the problem of the previous overestimation is nearly solved. Es- 
pecially, a drastic improvement is achieved for the S31 amplitude. The phase-shifts are better 
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Figure 7: Scattering amplitude for the Sn and S31 irN partial waves with the improved Cy. 
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Figure 8: Phase-shifts for Sn and S31 ttN scattering with the improved C^. 



reproduced now as shown in Fig|| These results show that the correction of eq. (|20"D is important 
and leads to improved results with respect to the plain use of the standard lowest order chiral 
Lagrangian. In the following sections, we employ this modified coefficient. The calculated inelas- 
ticities with this coefficient, are almost the same as before and show the lack of some important 
channels. We come back to this problem in the next section. 

3 ttN scattering in a 2-3-body coupled channel model 

In this section, we extend our model to include the irnN channels. The cross sections of the 
ttN — > tcttN scattering are known experimentally and they are sizeable compared with the two 
body cross sections. In this paper, we include only the transition potential between ttN and 
tttcN and disregard the coupling between the {KH, KA, rjn, nnN} states and the tttxN . This is 
a simplification forced by the ignorance of such couplings, but the larger mass of these states 
with respect to ttN, and the fact that we are talking about corrections, make this simplification 
justifiable. 

3.1 ttN — > ttttN process 

In the isospin formalism, the nN —>■ ttttN transition amplitudes can be classified by the total 
isospin I and the isospin of the final two pions I n7T , and the corresponding amplitudes are written 
as A 2 ii n7T . The amplitudes of the physical processes are expressed in terms of the following four 
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Figure 9: ttN — > ttttN process 



independent isospin amplitudes 



A n = an x\&- (K ~ k 2 )Xi , A 31 = a 3i x\°- (K - k' 2 )Xi 
Aw = a w Xf&- kxi , A 32 = a 32 Xf°' k Xi 



(22) 



where Xi an d Xf are spinors for the initial and final nucleon, and k , k[ and k! 2 are the momenta 
of the pions depicted in Fig.|9|. Among these amplitudes, the upper two amplitudes which have 
I nn = 1 correspond to transitions from tiN in an S-wave, while the lower two amplitudes correspond 
to the transition in P-wave. 

These amplitudes have been extracted by analyzing the experimental cross sections. In ref. 
, the data from threshold to 1470 MeV are analyzed. The reduced amplitudes a^ are assumed 



to be constant except for a 10 which has a resonance behavior due to a coupling to the Roper 
resonance, Pn(1440), 

«io(V5) = a' 10 M ~f. a r r ,2 (23) 

where \fs th ~ 1213 MeV is the threshold energy. The reduced amplitudes obtained are 
a u = 10.61 ± 0.62 [m^ 3 ] , a 31 = -6.02 ± 0.31 [m; 3 ] 



a 



i„= 6.63 ± 0.21 [m; 3 ] , a 32 = 2.75 ± 0.13 [m; 3 ] 



(24) 



with M = 1416 ± 14 MeV and T = 287 ± 43 MeV. In ref. El, the data close to threshold 



including newer data with respect to |23|, are analyzed. In this case the reduced amplitudes are 



assumed to depend linearly on the center of mass energy. The reduced amplitudes 

a u = 3.3±0.8 + (0.9±2.0)(^-^JK [m~ 3 ] (25) 
a 31 = -5.0 ±2.2 + (15.0 ± 4.2) (Vs-V%JMr [«V 3 ] (26) 
aio = 6.55 ±0.16 + (10.4 ±0.8)(Vs-Vs t h)/m* \ m ^\ ( 27 ) 
a 32 = 2.07 + 0.10 + (1.98 ±0.33) (Vi - y/s th )/m„ [m^ 3 ] (28) 

are obtained. 

We plot both empirical reduced amplitudes in Fig.|10|. One can see that the extracted ampli- 
tudes of S-wave irN are very different in both analyses, while those of P-wave irN agree quite well. 
These discrepancies reflect the difficulties to determine the amplitudes of the S-wave irN from 
(tt, 27r) cross section data, which looks quite natural since most of the processes are dominated 
by the transition from P-wave. Therefore, in this paper, we leave free the S-wave ttN amplitudes, 
an(\^s) and 031 (y/s), and determine them through the present study so that they are compatible 
with the data of ttN — > ttN elastic scattering. 
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Figure 10: Empirical amplitudes of ttN tcttN transition. The continuous lines and dotted lines 
correspond to the amplitudes of paper and respectively. The absolute value is plotted for 



3.2 The model including the 7tttN channel 

We consider the Bethe-Salpeter equation for the scattering matrix, eq. ([!]) with the eight coupled 
channels including ttttN, namely, ir~p , ir n , rjn , K + T,~ , K°Tf , K°A, 7r°7r~p and 7T + 7T~n. We 
do not include the 7r°7r°n channel because it does not couple to the S-wave nN state. 
The potentials of the ttN <-» ttttN transitions are written as 



^TT p,TT°1T p 



y/2 y/2 
-z-vn + —v-3,1 

3 



10 



Vr. 



1T~p , 7T~*~ 7T TL 



- 3 vn ~ ~ 3 v 31 

1 1 

--vn + -v 31 



X j f v-(k[-k' 2 )xi + 



10 

V2 Vh 



15 



^32 



x\°-kx 



XfV-(k[-k' 2 )xi ^-v^XfV-kxi 



y/2 y/2 



6 



-V31 



x j f °-(K-k' 2 )xi + 



--.10 



10 



15 



-V32 



(29) 
(30) 
(31) 

X^-kxi (32) 



in terms of the reduced potentials wh(a/s) , v 3 \(y/s), vio(y/s) and v 32 (y/s), which correspond to 
the reduced amplitudes an(yfs) , a 3 i(y/s) a io(v^) an d a 32(Vs) respectively, after solving the BS 
equation with all the channels. Note that, in our formalism, tt + corresponds to the — \I = 1,I Z = 1) 
state. The terms of v w (y/s) and v 32 (y/s) do not contribute in the present S-wave scattering. 



For the nivN state, we introduce the baryon-two meson propagator corresponding to Fig.|Tl 
' d 4 qi f d 4 q 2 



G{P) 



(2vr) 4 J (2tt 



:(qi - Q2)' 



2M, 



N 



1 



(P — gi — q 2 ) 2 - Mjf + it q\ - m\ + ie q 2 - ml + it 



(33) 



which includes the vertex structure for convenience. The integral in eq. (|33|) is strongly divergent. 
Its Imaginary part is finite and, neglecting the contribution of the negative energy of the baryon, 
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Figure 11: Propagator for the 7T7tA^ state. 
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Figure 12: Imaginary part of G(\/s). 



is given by, 
Im[G(v^)] 
where 



M N 



4(2vr) ; 



, du! jduj 2 \M 2 N + 2q 1 (co 1 ) 2 + 2g 2 (cu 2 ) 2 -{y/8-u x - u 2 ) 2 \ 9(1 - A 2 ) (34) 



.4 



(y/s-Uj- u 2 ) 2 - Ml - gi(^i) 2 - q 2 {uj 2 f 
2q l (ui)q 2 (u 2 ) 



(35) 



with qi(u>i) = Juj 2 — vn\. This imaginary part, reflecting the phase space for the intermediate 



baryon-two meson system, is shown in Fig.[12|. On the other hand, the real part of G(a/s), in the 
renormalized model, is not fixed. The situation is well illustrated in a dispersion relation 



dw- 



\m[G{w)} 



+ A + B(xfs - w ) 



(y/s-W ) 

7T A/i (W — W ) 6 (W - y/S) 

+C(yfs - w ) 2 + D{yfs - w ) 3 + E(yfs - wo) 4 + F(y/s - w ) 5 



where six arbitrary constants A ~ F are introduced because Im[G(\/s)] grows as 

lm[G(\/s)] j — oo for n < 5 



lim 

>oo 



-OO 





for 



n > 6 



(36) 



(37) 



which is easily proved by an actual calculation or dimensional considerations. Therefore, in this 
paper, we treat Re[G(\/i)] as a free function and look for reasonable values consistent with the 
experiment. 



3.3 Final results 

By varying the potential Vn(\/s) and Vsi(y/s), the subtraction constants cij(/i) and Re[G(\/s)], we 
try to reproduce the experimental elastic ttN T matrix and also the ttN — > ttttN cross sections. 
For the Vu(y/s) and Vsi(y/s) functions we employ a real polynomial function, and for Re[G(y / s)] 
we test several types of functions. 

In Fig.^ we show the propagators G(\/s) and G(\fs). The subtraction parameters ai(fi) for 
the meson-baryon propagators that we obtain in the new fit to the data, are 

fi = 1200 MeV, a nN {fi) = 2.0, a vN {fi) = 0.1, a KA {fi) = 1.5, a K ^{fi) = -2.8 (38) 

which are a little changed from the previous values omitting the tcttN channels. For the G(y/s) case 
the real part is taken zero. We find that the best results are obtained with functions compatible 
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Figure 13: Propagators of meson-baryon and tcttN systems using the parameters of eq 
solid and dashed lines are real and imaginary parts respectively 



The 



with zero. This result is similar to the one found in pHI where the real part of the three pion loop 
form decay was also found negligibly small. 

Fig.[F| shows the functions an( v / i) and a 3 i(y/s) determined in this study. In the upper two 
graphs, the potentials Vn(yfs) (f3i(\/s)) and the amplitudes an(y/s) (a,3i(\/s)), which come after 
unit arizat ion, are shown. The thin lines correspond to the following potentials which we finally 
use 

Vu(V^) = 4.0 + 1.0(^-1213)/^ [m; 3 ] (39) 
u 31 (Vs) = -5.60( v / s- 1470)/m^ - 1.05( v / s- 1470) 2 /^ + 1.77(v^- 1470)7^ (40) 
+0.66( v / s - 1470) 4 /mt - 0.17( v / s - 1470) 5 /^ - 0.07(v^ - 1470) 6 /™° [m. 



-31 



The continuous line and dashed line correspond to the real and imaginary part of an(y / s) (031 
respectively. One can see that the imaginary parts are almost negligible as assumed in the paper 
and P3|. In the lower two graphs, the real part of an(y/s) (a 3 i(y / s)) is plotted in comparison 
with the two empirical ones. The lower energy part of the calculated ou(i/s) agrees with that of 
the paper |24j but it is quite different from the one of the | 23fl . On the other hand our calculated 
03i(V^) amplitude is different from both [2^j and |^|. We should note that, while it is possible 
to reproduce the ttN — > 7r7rA^ cross sections with the three set of amplitudes, a simultaneous 
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Figure 14: The S-wave tcN <-> nnN amplitudes an(y/s) and 031(^/5). Upper figures in the panel: 
the thin lines are the potential %(\/S), the solid lines are Re[ay(\/s)] and the dashed lines are 
Im[ajj]. Lower figures in the panel: the calculated Re [a^] (solid) are compared with |ZB| (dashed) 
and @ (dotted). 



description of the irN — > ttttN cross sections and the ttN — > 7riV scattering date is not possible 
with the amplitude of p3| and |p4| . We shall elaborate further on this issue below. 



The T matrix elements obtained are shown in Fig. 15. One can see that the T matrix elements 



are well reproduced at energies below 1600 MeV in both the Sn and cases. The phase-shifts 
and inelasticities are shown in Fig.lTBl. The inclusion of the mrN channels improves slightly the 



phase shifts in Sn and S31. The important thing to note is that, as seen in Fig.|16|, the inelasticities 
at low energies in both Sn and S31 are well reproduced with the inclusion of the tttcN channels. 
The present ou(^/s) function is essentially the same as that of the paper [24| in the range of 
energies studied there. Should we take the amplitude of paper f23j, the inelasticities would be 
much overestimated at energies around 1400 MeV. The present a 3 i(y/s) amplitude is different 
from both empirical 031(^/5) determinations. It has a node at the energy 1470 MeV which is 
reflected in the inelasticities shown in the figure. The inelasticities are quite sensitive to the 
a 3 i( v /i) amplitude used. For example, the two empirical 031(1/5) correspond to the inelasticities 
shown in Fig.[T7] which differ appreciable from those obtained with the function determined here. 
We should note that for this test we have used a function i>3i(\/s) such that after unitarization 
they lead to scattering amplitudes identical to the empirical asi(y/s). The unitarization procedure 
modifies only little the potential as is visible in upper right box of Fig.[l4|for our case. The opposite 
sign of our 031 (-^/s) also provides the same inelasticities and hence is compatible with the S31 ttN 
scattering data. However, it leads to unacceptable results for the irN — * tcttN cross section. 

Fig.|l8] shows the ttN — > ttttN cross sections. They are calculated with the present S-wave 
amplitudes and the P-wave amplitudes of the paper |[23|| . The dashed lines correspond to the 
cross section when we drop the S-wave contributions. The two processes, ir~p — > 7r°7r°n and 
tt + P — > 7i + 7i + n, are purely P-wave ttN and have no S-wave 7iN contribution. In the 7r~p — > 7i + n^n 
reaction, the P-wave contribution is large and dominates the process. The effect of the S-wave 
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Figure 15: Scattering amplitude for the Sn and S31 irN partial waves with 7T7riV channels. 
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Figure 16: Phase-shifts and inelasticities of Sn and ixN scattering with 7T7riV channels. 
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Figure 17: Inelasticities of S31 ttN scattering which correspond to the empirical ttN — > ttttN 
amplitudes. The solids line in the left and right figure correspond to |23| and |2|| respectively. 
The dotted line is the experimental analysis 



is not visible in the figure. However, in the ir~p — > 7t°tt~p and n + p — > 7T°7T + p processes we can 
see that the P-wave contributions are small and do not explain the size of the data. As one can 
see, the present S-wave amplitude provides enough strength to account for the cross section data. 
In the n~p — > n°7i~p process, the S-wave amplitude is a mixture of an(y/s) and 031(^/5). In 
order to obtain the contribution of the S-wave shown in the figure, the amplitudes an(y/s) and 
03i(v^) should interfere constructively close to the ttttN threshold. The relative sign between 
axi(y/s) and a 3 i(-y/s) (t>ii(\/s) and 1*31 (\/s)) is determined by this condition. The value of an/a 31 
at threshold is about +1/2, which is in contrast to the one pion exchange prediction, —2 f31 |. 
However, more elaborate models for two pion production J32], [53|, [34], [35], |36|, |37] contain many 
more terms that change the value of this ratio. The amplitudes au(y/s) and a 31 (y/s) determined 
in this study, account for both the 7iN elastic scattering data and the tiN — > rntN cross section 
data simultaneously. 

In Fig.[L9] we also show the result for the Tt~p — > r/n cross section. We can see that the 



agreement is good up to 1550 MeV, above which higher partial waves, as demonstrated in p8 |, 
become relevant. 

We have also evaluated the scattering lengths in our approach which are listed in Table |2]. 
The thresholds for the K°A, K°TP, K+S-are = 1613, 1690, 1691 MeV respectively. These 
energies, particularly the two last ones, are already in the region of energies where the theory 
deviates from experiment in the phase shifts and inelasticities, therefore we should take these 
numbers only as indicative. On the other hand, the fitting to the data has been done around the 
iV*(1535) energy. Thus, our prediction for the r/n scattering length, a m = 0.26-H0.25 fm, should 
be rather accurate. This number is in agreement with the result quoted in 0, a vn = 0.20 + i0.26 
fm, although a bit more attractive. Still the real part is about a factor three smaller than in [21 
and a factor four smaller than in 
Re fa 



tjn J 



In spite of that, it was argued in |y that even the small value 
0.2 fm is not unrealistically small. This scattering length is important since it plays a 



crucial role in the possibility to have rj bound states in nuclei [40], 1411 . The scattering length for 



n n and n~p have not been imposed in the fit to the data, which concentrated in the iV*(1535) 
region, as already mentioned. In this sense, the agreement with the data about 450 MeV below 
that resonance should be considered an unexpected success. We obtain isospin 3/2 and isospin 
1/2 scattering lengths a 3 = — 0.0875 m" 1 , a\ = 0.1272 m" 1 , to be compared with the experimental 
numbers g2J, a e 3 xp = -0.0852±0.0027m~ 1 , af v = 0.1752±0.0041 m~ x . The agreement is good for 
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Figure 18: The cross section of ttN — > ttttN scattering. Data points are taken from the papers 



in reference of [24] 
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Figure 19: The cross section of ttN — > rjn scattering. 
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Table 2: The calculated meson-baryon scattering lengths in unit of [fm]. 





7r°n 


7T p 


rjn 




at [fm] 


-0.023 


0.080 + i 0.003 


0.264 + i 


0.245 




K°A 








^ [fm] 


-0.148 + i 0.165 


-0.205 + i 0.068 


-0.284 + i 


0.090 



the isospin 3/2 scattering length, but the isospin 1/2 one is about 25% smaller than experiment. 

To summarize the result of section 2 and 3, wee can see that the chiral unitary approach 
including the ttttN channels, is a very efficient tool to study the S-wave ttN scattering. It can 
reproduce, with few parameters, not only the isospin 1/2 part but also the isospin 3/2 part, 
which could not be obtained before [[| [19|, [20(1, when only meson-baryon channels were considered. 
Particularly, the irnN channel was found essential to reproduce the isospin 3/2 scattering. 



4 The iV*(1535) and A(1620) resonances 

As shown in Fig.^, the calculated T(Sn) has a resonance behavior around 1535 MeV indicating 
that the well known negative parity baryon iV*(1535) is generated in this approach. To find the 
pole corresponding to the resonance we extend our calculation of the T matrix to the complex Pq 
plane. We evaluate the T matrix elements by means of 



T(Po) 



V(Po)- 1 



G(P 



-i 



(41) 



and look for poles in the complex P° plane. In this plane the function G(Pq) has cuts on the 
real axis. For example, Fig.|2(] left shows Im[G 7r - p (Po)] m the physical sheet, namely the first 
Riemann sheet. We also plot to the right of Fig.^ the sheet connecting the first Riemann sheet 
for Im[Po] > with the second Riemann sheet for Im[Po] < 0, which is defined as 




2iIm[Gf(P )] 



at 
at 



Re[P ] 
Re[P ] 



< 
> 



(42) 



with y/s^i the i-th channel threshold energy . We need also to extrapolate G(y/s) to the complex 
plane. For this purpose we parameterize the result for Im[G(\/s)] above the imN threshold in the 
real axis as 

Im[G(Vi)] = -0.638( v / s - 1213)/m* + L124(Vi - 1213) 2 /^ - 0.882(v^ - 1213) 3 /^ (43) 

which allows for an analytical continuation for Re[Po] above that threshold. 

We search for poles of the isospin 1/2 T matrix elements on this sheet, and obtain a pole at 



P 1 

1 n 



1543 - % 46 MeV 



(44) 



The structure of the T matrix guarantees that we find the pole in all the elements which have an 
isospin 1/2 component. For example, Fig.^1] left shows |T(5ii)|, where we can see a pole clearly. 
This result tells us that the decay width of the resonance is about 93 MeV. This value is smaller 
than the PDG estimation 100 ~ 250 MeV g|, but agrees with the new data from BEPC, 95 ± 15 
MeV p|. 
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Figure 20: The imaginary part of the propagator for the it p system. 




Figure 21: |T(Sii)| (left) and |T(S3i)| (right) in the complex energy plane. The poles of iV*(1535) 
and A(1620) are seen at the energies 1543 — i46 MeV and 1625 — z215 MeV respectively. 



On the other hand, in the elements of pure isospin 3/2 we find a pole around 1625 — z215 MeV, 
which is reminiscent of the A(1620) resonance, although with a larger width than the nominal 



one of the PDG of about 120-160 MeV. In Fig. 21 right we show |T(S f 3i)|, where a pole is clearly 



visible. We should consider the position of the pole only as qualitative in our approach since 
around yfs = 1600 MeV is where we start having noticeable discrepancies with the data and 
furthermore we had to go far in the complex plane where the analytic extrapolation of the mrN 
propagator becomes less accurate. This resonance is responsible for the change of curvature in 



ImT(5 , 3i) shown in Fig. 15. It is interesting to note that the introduction of the nirN channels is 
what has induced the appearance of this resonance, which does not show in our approach when 
the ttttN channels are not considered. In this respect it is interesting to note that the A(1620) 
resonance couples mostly to the irnN channel [£| . 



For Pq near the pole the T matrix elements are approximated by 
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Table 3: Coupling constants and decay widths of iV*(1535). 











K°A 


71 p 


7r°n 


rjn 


7T 7T p 


7T + 7T U 




\9i\ 


2.12 


1.50 


0.92 


0.56 


0.39 


1.84 


0.57 m~ 2 


0.40 m~ 2 


r 


\ [MeV] 








14.1 


7.0 


65.7 


4.6 


2.4 


B. 


ratio [%) 








15.0 


7.4 


70.1 


4.9 


2.5 




9i 


-2.36 


1.67 


-1.28 


-0.57 


0.39 


1.77 


-0.61 m~ 2 


-0.43 m~ 2 


I 


\ [MeV] 








15.0 


7.1 


60.8 


5.4 


2.6 


B. 


ratio [%] 








16.5 


7.8 


66.9 


5.9 


2.9 



with gi the couplings of the resonance to the i-th channel. We obtain the size of the couplings by 
evaluating the residues of the diagonal elements Ta(P ). The values obtained are listed in Table 
|3|. We find that the resonance iV*(1535) couples strongly to the KY> and rjn channels strongly 
with a couplings four times bigger than for the ttN channels. The couplings to the KA are also 
large compared to the tiN channels. In short, we get 

\9ttn\ < \9ka\ < \g V n\ ~ \9kt\ ■ (46) 

Using these couplings we calculate the partial decay widths of the open channels by means of 

T^^lmlG^M^M 2 (47) 

where we take Mn* = RefP^] = 1543 MeV. The partial decay rates and the branching ratios 
obtained are also listed in Table |3|. The calculated branching ratios of ttN, t]N and ttttN decay 
modes, are 22%, 70% and 7% respectively. The fraction of rjN mode is large, which is know to be 
characteristic of this resonance although our fraction is bigger than the PDG estimation 30 ~ 55 
%. Our ttN fraction is smaller than the PDG estimation, 35 ~ 55 %, while the ttttN fraction is 
compatible with the PDG estimation, < 10 %. 

The lower part of Table [3] shows the same quantities obtained from a Breit-Wigner fit of the 
real energy scattering amplitudes T^y'i). We fit them by the Breit-Wigner form together with a 
background 

TijiVZ) = Mj ^_ 1543 + i46 + % + hAV^ ~ 1450) (48) 

at 1450 < ^/s < 1650 MeV. The unknown parameters ft x gj, and 6^ are determined by the 
method of least squares. We obtain the values of gi from the g vn x g^ corresponding to the rjn final 
state amplitudes. Their absolute value agree fairly well with the values obtained from the pole 
residues and gives us confidence about our numerical evaluation of the couplings. For instance, 
the ttN, rjn and ttttN branching ratios are now 24 %, 67 % and 9 % respectively. 

This latter analysis allows us to determine the sign of the couplings. The signs given are 
relative to that of g m . It is instructive to decompose the resonance in the SU(3) representations. 
In fact, our result, ignoring the coupling to the titiN channel, leads to 

g 8 = -2.52 , g w = 2.62 , g l0 = 0.43 , g 27 = -0.47 (49) 

and tells us that the iV*(1535) resonance is almost an equal weight mixture of the R-parity even 
SU(3) octet 8 and the R-parity odd SU(3) octet 8'. It would be interesting to compare these 
results with the results from other models or lattice QCD simulations. 
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5 Conclusion 



We have studied the S-wave ttN scattering, together with that of coupled channels, in a chiral 
unitary model in the region of center of mass energies from threshold to 1600 MeV. We calculated 
the T matrix using the Bethe-Salpeter equation in the eight coupled channels including six meson- 
baryon channels and two nirN channels. We took the transition potentials between the meson- 
baryon systems from the lowest order chiral Lagrangian and improved them taking into account the 
vector meson dominance hypothesis. Then we introduced the appropriate 7rN <-> ttttN transition 
potentials which influence both the elastic scattering and the pion production processes. In the 
present model the renormalization due to higher order contributions is included by means of 
subtraction constants in the real part of the propagators of the two or three-body systems, which 
are taken as free parameters and determined through comparison with the T matrix of the data 
analysis. The imaginary part of the meson-baryon or ttttN propagators is fixed and ensures 
unitarity in the S matrix. A realistic T matrix is obtained with a few free parameters for energies 
up to 1600 MeV. The phase-shifts and the inelasticities are well reproduced in both isospin 1/2 
and 3/2. We find that the correction of the chiral coefficient and the 7nrN channels are important 
to obtain an accurate T matrix, especially in isospin 3/2. Our analysis allowed us to determine 
the S-wave amplitudes for ttN — > rnrN and we found that the isospin 3/2 nN <-> rniN amplitude 
is different from the two previous empirical ones. 

The resonance AT* (1535) is generated dynamically and qualifies as a quasi bound state of meson 
and baryon. The corresponding pole is seen in the T matrix on the complex plane. We calculate 
the total and partial decay width of the resonance. The total width obtained, about 80 MeV, is 
smaller than the PDG estimation, but agrees with the new data from BEPC. Also the large r\n 
branching ratio observed in the data is reproduced. 

The present study has served to show the potential of the chiral unitary approach extending 
the predictions to higher energies than it would be possible with the use of xPT . Yet, we also 
saw that improvements in the basic information of the lowest order chiral Lagrangians to account 
for phenomenology of VMD are welcome. On the other hand we found mandatory the inclusion 
of the 7T7iN channels in order to find an accurate reproduction of the data, particularly those 
in the isospin 3/2 sector. We also found that the introduction of these channels, forcing them 
to reproduce the inelasticities and other data, has as an indirect consequence that the A(1620) 
resonance appears then as a pole in the complex plane indicating a large mixing of this resonance 
with 7mN states. This interpretation would be consistent with the large experimental coupling of 
this resonance to the nnN channel. 
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